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We study the critical point for the emergence of coherence in a harmonically trapped two-dimensional Bose 
gas with tuneable interactions. Over a wide range of interaction strengths we find excellent agreement with 
the classical-field predictions for the critical point of the Berezinskii-Kosterlitz-Thouless (BKT) superfluid tran¬ 
sition. This allows us to quantitatively show, without any free parameters, that the interaction-driven BKT 
transition smoothly converges onto the purely quantum-statistical Bose-Einstein condensation (BEC) transition 
in the limit of vanishing interactions. 
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Reducing the dimensionality of a physical system increases 
the importance of thermal fluctuations and can profoundly 
affect the type of order that the system can display at low 
temperatures GM- In a uniform two-dimensional (2d) 
Bose gas, the true long-range-order (LRO) associated with 
Bose-Einstein condensation (BEC) is precluded at any non¬ 
zero temperature by the Mermin-Wagner theorem. However, 
an interacting 2d Bose gas still undergoes the Berezinskii- 
Kosterlitz-Thouless (BKT) transition to a superfluid state at 
a critical phase-space density Ubkt 0161 . Quantitatively 
predicting Ubkt requires an accurate description of the non- 
perturbative behaviour of interacting bosons in the fluctuation 
region near the critical point. Classical-field simulations Q 
model this behaviour by a turbulent matter-wave held ii, and 
predict Ubkt = In (380/p), where p is a dimensionless mea¬ 
sure of the interaction strength 0. This result makes it man¬ 
ifest that the transition is interaction driven; the critical tem¬ 
perature, Tc oc n/DsKT, where n is the gas density, vanishes 
in the non-interacting limit, g —t 0, for any non-infinite n. 
While for phase-space density D > Ubkt true LRO is still 
absent, the first order correlation function, pi(r), decays only 
algebraically at large distance, in contrast to the exponential 
decay in the normal degenerate gas. Such extended coherence 
is sufficient for superfluidity Qol, and in practice offers a sig¬ 
nature of the phase transition lHIHIl- 

In contrast to the infinite uniform system, in a 2d harmonic 
trap, pertinent to most ultracold-atom experiments on BKT 
physics iri5H29l . the modified density of states allows for a 
BEC transition to occur in the ideal gas (p = 0) Eolllll (see 
also 1321). In an isotropic trap of frequency oJr, it should occur 
at a critical atom number BOl [311 
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This ideal-gas BEC transition is similar to the familiar con¬ 
densation in 3d; it is purely quantum-statistical and follows 
from Einstein’s standard argument of the saturation of the 
excited states. However, the importance of dimensionality 
emerges if interactions are introduced. In 3d, ideal-gas BEC 
occurs when in the trap centre D « 2.612; weak interactions 


slightly shift the critical atom number Il33l - l35l . but do not al¬ 
ter the BEC-like nature of the transition. In 2d, while iV° is 
finite, the phase-space density required in the trap centre for 
the excited states to saturate is infinite, just as in a uniform 
system na. Eor any p > 0 this is unattainable and the ex¬ 
cited states can accommodate any number of particles 1^ . 
The BEC transition is thus suppressed and one expects it to be 
replaced by the BKT transition with non-infinite Ubkt m. 

These arguments suggest that the two conceptually very dif¬ 
ferent phase transitions, the interaction-driven BKT and the 
saturation-driven BEC, are in fact continuously connected as 
g —?■ 0 IIToIIJTI . The harmonic trapping potential offers the 
opportunity to experimentally observe this unification of BKT 
and BEC physics. While in an infinite uniform gas no tran¬ 
sition occurs for p = 0, in a harmonic trap a transition al¬ 
ways occurs at a finite Nc and always results in significantly 
extended coherence of the gas. The nature of this transition 
is quantitatively encoded in the value of Nc- This picture 
is supported by the calculations of the critical atom number 
|[T0]|38l, based on the classical-field simulations Q[33, 
which suggest that smoothly connects to TV® in the limit 
p = 0. In a finite-size system the change from the BKT to the 
BEC transition is a crossover that spans a nonzero range of p 
values, but for realistic experimental parameters the width of 
this crossover region is very small, p < 10’ jUSlIlOl. 

Various signatures of a phase transition, including emer¬ 
gence of extended coherence lITSl [T 7 UT 9 I I2TII and superfluid¬ 
ity EM . have been observed in trapped 2d gases with different 
specific p values. On the other hand, systematic studies with 
tuneable p have focussed on in-situ measurements of the equa¬ 
tion of state ||22l|27l, which do not directly reveal any striking 
signatures of the infinite-order BKT transition. 

In this Letter, we systematically study the critical point 
for the emergence of extended coherence in a harmonically 
trapped 2d Bose gas over a wide range of interaction strengths, 
0.05 < p < 0.5. We show, without any free parameters, that 
Nc generally agrees very well with the beyond-mean-fleld cal¬ 
culation of li3^ . and converges onto iV® of Eq. ([^ as 
p —0. The critical chemical potential, jjtc, which directly re¬ 
veals uniform-system conditions for a phase transition to oc- 
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cur in the trap centre, also agrees with the BKT theory and 
converges onto the BEC value, pc = 0, for g —i' 0. Our mea¬ 
surements also reiterate the importance of the suppression of 
density fluctuations in the normal state near the BKT critical 
point, previously observed in lfT8l[T9ll2TH23]l . 

The experiment was carried out using a gas, in the 
apparatus described in HD- For 2d trapping, the tight axial 
(vertical) confinement is provided by two repulsive “blades” 
of blue-detuned light, formed by passing a 532-nm gaus- 
sian beam through a O-tt phase-plate MM, while a red- 
detuned 1064-nm dipole trap provides the in-plane (horizon¬ 
tal) confinement. The radial and axial trapping frequencies 
are {uJr,LOz) « 27r x (38,4100) Hz. For all our measurements 
T G [140 nK, 190 nK] and fi/ks < 100 nK, resulting in a 
small (< 30%) occupation of the excited axial states. The in¬ 
teraction strength, g = lfT4l . where a is the s-wave 

scattering length and tz, = [rnujz), is controlled via a 
Feshbach resonance centred at 402.5 G umia. 

To characterise long-range coherence of a gas we study its 
(in-plane) momentum distribution n(fc) lfT9l . A change in 
the functional form of gi (r) leads to a dramatic change in its 
values at distances much larger than the thermal wavelength 
A = h/y/2TTmkBT iflTll . and an increase of coherence over 
some large distance L manifests itself in enhanced population 
of the low-momentum states k ^ 21: jL. Thus, unlike the in¬ 
trap density distribution, which varies very smoothly through 
the BKT critical point ifTSl [191 l22l |2^ . n(k') can provide a 
dramatic signature of the phase transition US). 

As illustrated in Fig. [T] to identify the critical point for a 
given g, we start with a highly coherent 2d gas and mea¬ 
sure n(fc) after holding the cloud in the trap for a variable 
time t. During the hold time, the atom number N slowly 
decays through various inelastic processes m, while the 
elastic-collision rate (« 0.2 Ng'^ s“^) remains sufficiently 
high to ensure that the gas is in quasi-static equilibrium. To 
measure n{k), we employ the “momentum focussing” tech¬ 
nique imiiiiiiiiiii. We turn off just the tight 2 ; confine¬ 
ment, so the rapid vertical expansion (predominantly driven 
by the zero-point motion along z) removes all the interaction 
energy on a timescale 1/uJz ^ l/uJr- The subsequent hori¬ 
zontal ideal-gas evolution in the remaining in-plane harmonic 
potential reveals n{k) as the spatial distribution after a quarter 
of the trap period. We probe this distribution by absorption 
imaging along 2 : [see Fig.[2a)]. 

Our fc-space imaging resolution, Afc « 0.4 /rm“^, sets the 
largest distance over which we can probe coherence to L = 
27r/Afc « 15 pm, which is much larger than A « 0.7 pm. To 
probe coherence on this lengthscale, we simply monitor the 
peak value of the momentum distribution, Pq, without mak¬ 
ing any theoretical assumptions about the exact shape of n{k) 
at low k. To get the corresponding atom number A we do a 
simple summation over the image. Importantly, we eliminate 
the systematic error due to the uncertainty in the absorption¬ 
imaging cross section by independently calibrating our imag¬ 
ing system through measurements of the BEC critical point in 
a 3d gas HtI . 
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FIG. 1. (color online) Determination of the critical point for the on¬ 
set of coherence, for g = 0.28 and T ~ 140 nK. (a) Evolution of the 
momentum distribution, n{k), with the hold time t (see text). Ex¬ 
tended coherence is revealed as a sharp peak in n{k). Each image 
is an average of three experimental realisations, (b) Evolution of the 
momentum-distribution peak, Pq, and the smoothly decaying total 
atom number N. We associate the threshold-like behaviour of Pq 
with the critical time and deduce the corresponding Nc- The solid 
line is a heuristic piecewise fit function used to determine (m. 


In Fig. [^b) we show a typical evolution of Pq and N (here 
g = 0.28). While N decays smoothly, Pq shows two distinct 
regimes, which allows us to identify the critical hold time tc 
and the corresponding Nc- We note that even for N signif¬ 
icantly below Nc the peak of n{k) rises above a Gaussian 
fitted to the wings of the distribution, indicating some coher¬ 
ence on a lengthscale > A imiii]. The smooth evolution 
of such non-Gaussian “peakiness” of n{k) does not reveal a 
phase transition liSTIl . and only Pq corresponding to L ^ A 
shows a clear change in behaviour at a well-defined Nc ED. 
Our large L is still small compared to the thermal diameter 
of the cloud, 2y/kBT / (muj^) « 50 fim, so the observed Nc 
is closely linked to the occurrence of a phase transition in the 
centre of the trap 1471 . 

For comparisons with theory, we also fit p and T to each 
n{k) image. Unlike in 3d, in 2d interactions affect n{k) ap¬ 
preciably even in the normal state, and near the critical point 
it is in general insufficient to treat them at a mean-field (MF) 
level. However, beyond-MF correlations primarily affect the 
highly-populated low-A: states l38l . We restrict our fits to the 
high-fc wings of the distribution (U^k^ > 2gmkBT), where 
we expect the beyond-MF effects to be small, and still care¬ 
fully include the effects of interactions at a MF level l47l . 
Following E3, we also account for the residual thermal oc- 
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FIG. 2. (color online) Critical atom number as a function of the inter¬ 
action strength g. All numbers are scaled to the ideal-gas BEC crit¬ 
ical number A(?, defined in Eq. Solid line is the classical-field 
BKT prediction of Eq. l|^, without any free parameters. Dashed line 
is an approximation which neglects suppression of density fluctua¬ 
tions in the normal state. The star (*) denotes the critical point for 
BEC, which only occurs in the ideal-gas limit. The shaded region, 
g < 0.06, indicates the regime in which our measurements stop be¬ 
ing reliable (see text). Error bars are statistical. 


cupation of the axial excited states and the interaction-induced 
deformation of the axial eigenstates. 

In Fig. 1^ we summarise our measurements of the criti¬ 
cal atom number for a wide range of interaction strengths. 
To compare our data with the strictly-2d theoretical calcula¬ 
tions, we correct the observed “raw” Nc by subtracting the 
calculated population of the excited axial states Il47l . We 
scale this corrected critical number, Nc, to the BEC critical 
atom number, iV° of Eq. ([^, and plot it versus g. Our Ak- 
limited value of L imposes a lower bound on g for which 
we can reliably identify the critical point. In the absence of 
any phase transition, in a weakly-interacting degenerate gas 
gi{r) ~ exp(—r/£o), with fg = Aexp(I?/2)/'\/4^ ifTdl . 
We thus do not expect our measurements to reliably iden¬ 
tify Nc if io > L for some D < Hbkt- This occurs for 
g < 380 A^/(TttL^) « 0.06, indicated by the shaded area in 
Fig. I Our measurements also stop being reliable for 5 > 0.5; 
in that regime our ME temperature fits are restricted to very 
high k values, which are affected by the anharmonicity of the 
optical trap. The error bars in Eig. |^are statistical, while the 
systematic uncertainty in iVc/At° is < 0.2 1471 . 

Without any free parameters, we find generally excellent 
agreement with the prediction of Ref. ll^ : 
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which is based on fixing the phase-space density in the trap 
centre to D^kt and integrating a uniform-system equation of 
state over the trap, using classical-field results of Ref. IWI . 

The agreement with Eq. Q over a very broad range of in¬ 
teraction strengths, and the proximity of our lowest reliable g 


EIG. 3. (color online) Critical chemical potential as a function of the 
interaction strength g. Solid line is the classical-field BKT predic¬ 
tion. Dashed line is derived assuming a fully fluctuating Bose gas, 
while dotted line corresponds to a complete suppression of density 
fluctuations. Error bars are statistical. 


values to zero, allow us to conclude that the critical atom num¬ 
ber, without any free parameters, indeed smoothly converges 
onto the BEC result of Eq. Q- 

It is instructive to also compare our data with the approxi¬ 
mation N^^'^/Nc = 1 + 3 piTgj;j/ 7 r^ lll 0 ][T 2 l, shown by the 
dashed line in Eig. Here, the critical phase-space density is 
again set to Z^bkt, but the suppression of bosonic fluctuations 
in the normal state is neglected, i.e. the density profile is cal¬ 
culated using ME theory with an interaction potential 2gn{r), 
where g = /m)g. Our data strongly exclude this result, 

confirming the importance of the suppression of density fluc¬ 
tuations near the critical point even for our lowest g values. 

Eor a more direct comparison with the uniform-system the¬ 
ory, we also consider the critical chemical potential for the 
onset of coherence. Like Nc in Eig. [T] g,c is experimen¬ 
tally defined via the critical hold time tc- The classical- 
field simulations Q predict Hbkt to be reached for = 
k-NT {(j/tt) In (13.2/p), which reduces to the BEC prediction, 
fic = 0 , for g = 0. 

In Eig. we plot pc = versus g, and again ob¬ 

serve generally good agreement with the classical-field pre¬ 
diction (solid line), all the way down to p « 0.06, i.e. very 
close to the expected BEC limit. The small systematic dif¬ 
ference between the data and the theory is comparable to our 
systematic uncertainty in jlc, of ~ 0.05 El. 

We also compare our data with two intuitive approxima¬ 
tions to fic- We consider interaction potentials 7 pn with 
7 = 2, corresponding to a fully fluctuating Bose gas, and 
7 = 1 , corresponding to a complete suppression of den¬ 
sity fluctuations. In both of these extremes one can analyti¬ 
cally write Dj (/t) = — In [1 — exp (p — 7 pn/ (k^T))] ifT^ . 
Defining /(/ so that D,y{fi2) = D^kt we obtain the dashed 
(7 = 2) and dotted (7 = 1) lines in Eig. Generally 7=1 
provides a better approximation, highlighting how strong the 
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suppression of density fluctuations in the normal state is. 

Finally, we note that in previous experiments ll22l l27ll . on 
the in-trap equation of state, jlc was deduced by defining it so 
as to satisfy the theoretical expectation 17] [53 that the phase- 
space density should be a universal function of (/i — fj,c)/g- 
Our measurements of flc defined through the emergence of 
extended coherence show that the two definitions indeed lead 
to very similar values. 

In conclusion, by studying the critical point for the emer¬ 
gence of extended coherence in a harmonically trapped 2d 
gas with tuneable interactions, we have quantitatively con¬ 
firmed the predictions of classical-field theory and observed 
the expected unification of BKT and BEC transitions in the 
limit of vanishing interactions. The in-plane harmonic po¬ 
tential enables this observation. However, to quantitatively 
study the exact functional form of the slowly decaying cor¬ 
relations in a BKT superfluid, in the future it would be very 
interesting to study coherence of a tuneable 2d gas in a uni¬ 
form potential 1152145^ . Just below Tc this should reveal an 
interaction-strength-independent algebraic decay of the first- 
order correlation function, corresponding to a universal jump 
in the superfluid density li55l . 
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